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We show how the t+t~ spin correlations can be used to improve the 
recognition of the parent boson spin, and hence to identify scalar boson 
TJ" T^T^ events from the vector boson t^t^ background in 

high energy accelerator experiments. 

1. Introduction 

The most unsatisfactory feature of the Standard Model is our lack of 
knowledge of the actual mechanism that breaks the electroweak gauge sym- 
metry and generates the particle masses. In the Standard Model the break- 
ing occurs by a complex Higgs boson doublet. Three components of this 
doublet become the longitudinal polarisation states of the massive vector 
gauge bosons (VF^, Z), while the remaining component manifests itself as 
a neutral massive scalar Higgs particle. On the other hand the minimal 
supersymetric model contains two Higgs boson doublets. Three fields are 
taken by the vector bosons and remaining five become physical particles: 
a pair of charged boson H^, two scalars hP, and one pseudoscalar A^. 
Many more complicated mass generation models have been proposed, but 
their common feature is that the couplings of the Higgs bosons to fermions 
increase with the fermion mass. Thus for leptons the Ht^t~^ coupling dom- 
inates over the HiJ,~^fj,~, He~^e~ couplings. The experimental observation 
which breaks e, fi, r, universality by favouring t~^t~ events could be an 
indication of the presence of a Higgs scalar. Whenever such a departure 
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from lepton universality is observed a simple helicity correlation test can be 
performed which will clearly indicate that the r+r~ pairs have scalar boson 
origin as opposed to arising from vector boson decay. We will show that 
correlation of polarisations of the outgoing r leptons are very different for 
— ^ T~^T~ decays from that for Z/7* — > t~^t^ decays. The distinction 

arises because the vector bosons Z/j* decay into either + + or whereas 

the scalar Higgs bosons H into H — or — h, where +, — denotes the r pairs 
spin configurations. 

2. Universal interface for TAUOLA package 

We present the algorithm for the interfacing the r lepton decay pack- 
age TAUOLA - a library of Monte Carlo programs to simulate decays of r 
leptons [1-3] with "any" production generator to include spin effects in the 
elementary Z/'j* t^t~ process [4]. The approximate spin correlation are 
calculated from the information stored in the HEPEVT common block [5] filled 
by "any" r production program as described in Ref. [6]. As a demonstration 
example the interface is combined with the JETSET generator, however it 
should work in the same manner with the PYTHIA [7], HERWIG or ISA JET 
generators as well. In fact, such an interface can be considered as a sep- 
arate software project, to some degree independent both from the specific 
problem of r production and its decay. 

The aim of this interface is not to replace the matrix element calculations, 
but rather to provide a method of calculating/estimating spin effects in 
cases when spin effects would not be taken care of, at all. Such an approach 
is limited for the treatment of longitudinal spin degrees only and to the 
case of particle production and decay in the ultra-relativistic limit. The ap- 
proximation consists of reconstructing information of the elementary 2^2 
body process // (Z/7*) t^t~, hurried inside multi-body production 
process such as for example // — gZ, ff — 7Z, fg — fZ, — > fZ 
etc. The additional particles are grouped (summed) into effective quarks 
and leptons to minimise their virtualities. Such an approach is internally 
consistent in the case of photon or gluon emission within the leading log ap- 
proximation. The principle of calculating kinematic variables is simple. The 
4-momenta of the 2 — > 2 body process have to be found. The 4-momenta 
of the outcoming r's are at present used directly ^. Initial state momenta 
are constructed from the incoming and outcoming momenta of particles ac- 
companying production of the state Longitudinal polarisation of 
T leptons Vr depends on the spin quantum number of the r mother. It is 
randomly generated as specified in Table 1. The probability Pz used in the 

^ This part of the algorithm will be improved in a near future. 

^ The Z/f* state does not need to be explicitly coded in the HEPEVT common block. 
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Tabic 1. Probability for the configurations of the longitudinal polarisation of the 
pair of T leptons from different origins. 



Origin Vt-+ Vt- Probability 



Neutral Higgs bosons: h^,H^,A^ 




= +1 




= -1 


0.5 






= -1 


Vr- 


= +1 


0.5 


Charged Higgs boson: or H~ 




= +1 


Vr- 


= +1 


1.0 


Charged vector boson: or W~ 




= -1 


K- 


= -1 


1.0 


Neutral vector boson: Z/j* 


Vr+ 


= +1 


K- 


= +1 


Pz 




Vr+ 


= -1 




= -1 


l-Pz 


Other 


Vr+ 


= +1 


K- 


= +1 


0.5 






= -1 


Vr- 


= -1 


0.5 



generation, is calculated directly from the squares of the matrix elements of 
the Born-level 2 — > 2 process // — > r~r+: 

\M\)j^,-,^ (+,+) + \M\)f^,-,^ (-,-) ' 

where f = e, /j:, u, d, c, s, b. It can be also expressed (following conventions of 
Ref. [8]), with help of the couplings of fermions to the 7 ( and Z) bosons. 



f^(s,cose;l) + fe(s,cose;-l 
(s, cos cos 9Fi(s 



' ^ dcTBorn ( „ fl. 1 _1_ 5[£Bo-" ' ^ - X V ^ 

t^O'Born ^ ^ „„„ fl. „^ _ ^^ i „2 

dcos Q 

-p[(l + cos^ B)F2{s) + 2 cos ^F3(s)] . (3) 

with the four form-factors 



2 

i^o(s) = ^(9/9r + 2Rex(s)g/g^i;/t;^ + \x{s)f{v) + a))(vl + a^)), 

+ |x(s)|' 2t;/a/ 2^;^a^), (4) 



vra 



2 



^2(5) = ^( 2Rex(s)g/gr'y/ar + |x(s)r(^^/ + a/) 



4 



and 

'^''"^ ~ s - M| + isTz/Mz 



The Qf, Vf, af, q-,-, Vt, a,- are the charges and Z couphng constans of the 
fermions and r respectively. 




Fig. 1. Single n energy spectrum in the case of t produced from H (left-hand side) 
or Z (right-hand side), ^/s = mn or ^/s = mz respectively. 



3. Spin sensitive observables 

The polarisation of the r lepton can be exploited to identify a neutral 
Higgs bosons via the decay t'^t~ . Any experimental observation 

which breaks e, jU, r universality - the equality of e, ^ and r couplings to 
the gauge bosons - by favouring t~^t~ events could be an indication of the 
presence of a neutral Higgs scalar. The couplings of the "Higgs" parti- 
cles to fermions increase with the mass of the fermion, thus the r couples 
preferentially in comparison with either ji or e. Here the background is 
Z/7* — > r+T~ decays. Whenever such a departure from lepton universality 
is observed there exists a simple polarisation correlation test which, if used 
will help to indicate the presence of r pairs of Higgs boson origin among 
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(arbitrary units) 



(arbitrary units) 



Fast- Fast 
Slow-Slow 



Fast-Slow 
Slow-Fast 



Fast- Fast 
Slow-Slow 



Fast-Slow 
Slow-Fast 



Fig. 2. TT+TT energy- energy correlations in the case of t produced from H (left-hand 
side) or Z (right-hand side), ^/s = uih or ^/s = mz respectively. 



the background Z/j* — *■ r+r~ events. As we can see from Table 1, the r 
pairs arc produced with the well defined spin configurations: H — h or — — 

for vector bosons H or — |- for neutral Higgs boson. Thus a polarisation 

correlation test can be performed using the energy distributions of the fi- 
nal decay products which are sensitive to the admixture of the to Z/ 7* 
parentage of the r pairs. The polarisation correlation can be studied using 
the various r decay modes. Let us concentrate on the case of r decays to 
TTu, most sensitive to the spin correlations. The leptonic decay mode i.e. 
T~ e~Uci'-r- T~ ^ fJ'~'^ii'^T arc not very sensitive to the polarisation cor- 
relation because of the two missing neutrinos in this decay. In the case of 
other hadronic r decay, e.g. t~ — >■ i'tP~ i'tT^~'k^, t~ — ^ UtUi — ^ z/T-7r~27r°, 
T~ i't-K*~ Ut-K~tt^ etc. the reconstruction of vr^ is also necessary. For 
the production of the r lepton pairs Monte Carlo program PYTHIA was used, 
and for the decay Monte Carlo program TAUOLA, and this spin interface [4]. 
It was assured, that the invariant mass of the pair of two incoming quarks 
was ^/s = mz = mn- Energies arc defined in the t~^t~ pair rest-frame. 
With the help of variables z± = 2E^±/y/s, the spin effects are visualized. 
Fig. 1 shows the slope of tt energy spectrum (in the case of Z) due to the 
r polarisation. The slope of the distribution is simply proportional to the 
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polarisation. 

-^~1 + :P^2(z±-0.5) (6) 

dz± 

In the case of the plot on the left-hand side the spectrum is flat, as would be 
in the case of scalar neutral Higgs boson where there is no polarisation. As 
we can see in the Fig. 2 (vr+vr" energy-energy correlations) in Z/^* t^t~ 
decays a Fast (a Slow) tt^ is most likely to be associated with a a Fast (a 
Slow) TT^ whereas the opposite is favoured for — > t'^t~ decays. There- 
fore, a excess arising from Higgs boson decay can be recognised in the vr+vr" 
mode as a Fast tt^ with a Slow tt^. The quantity which can be measured 
experimentally is the invariant mass distribution. Fig. 3 shows tt'^tt" invari- 
ant mass distribution for the Higgs boson and Z cases. Continuous line - 
with spin effects included, dotted line - with spin effects switched off. Left- 
hand side plot corresponds to the Higgs boson case, right-hand side to the 
Z. In the case of Higgs boson, the mass distribution is peaked centrally, 
whereas in the case of Z/j* shoulders of the distributions are more pro- 
founded. In the Z/'j* case the Fast-Fast and Slow-Slow configurations are 
localised mostly at the shoulders of the tt'^tt'" invariant mass distributions, 
while for the the Higgs boson case the Fast-Slow configurations are localised 
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Fig. 4. The tt+tt" invariant mass distribution for uu Z/j* (left-hand plot) and 
dd — > Zf-y* (right-hand plot) produced with cms energy of 300 GeV. Continuous 
line with spin effects included, dotted line with spin effects switched off. 

in the centre of the distributions. If all polarisation effects are switched off 
(dashed lines) the distributions in the two cases are identical. This observ- 
able, i.e. a well defined distribution of invariant mass built from the visible 
decay products of the r's, can be helpful in separating Higgs boson signal 
from Z/^'* background. The same distribution have also been studied for 
the off-peak production of Z/^'* , i.e. for the larger cms energies. In these 
cases the average polarisation is large and negative, also distinct for the uu 
and dd annihilations. This may open a way for measuring the flavour of 
the quarks leading to r pair production. As illustrated in Fig. 4, the effect 
on the tt'^tt" invariant mass distribution is noticeable. The shape of the 
distribution might give the insight to the structure functions of colliding 
protons. 

4. Case of the Higgs boson signatures at LHC 

The r leptons are considered as a very promising signature for the 
searches of the Higgs bosons in the Minimal Supersymmetric Standard 
Model (MSSM) at LHC collider [9,10]. The neutral Higgs bosons H and 
A decay into t'^t~ pair, are enhanced for the large values of tan/3 (tan/3 
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Fig. 5. The tt+tt" invariant mass distribution after basic selection. On the left for 
Higgs boson, on the right for Z/'y*. Continuous line with spin effects included, 
dotted line with spin effects switched off. The Higgs boson mass was assumed to be 
300 GeV (see text). 



denotes the ratio of the vacuum expectation values of the Higgs doublets 
in the MSSM model), with the branching ratio of about 10% for most of 
the range of the interesting Higgs boson mass values ( 150-1000 GeV ). The 
irreducible background to this process is a Z/7* tt decay and the re- 
ducible backgrounds is the QCD jets. In this study the r identification is 
based on the presence of a single hard isolated charged hadron in the jet 
using tracker information. Two hard tracks from t~ and r"*" in the signal 
events have an opposite sign while no strong charged correlation is expected 
for the QCD jets. Fig. 5 shows the effect of tttt invariant mass distribution. 
Events were generated with the Monte Carlo PYTHIA 5.7 [11], the Higgs 
boson mass of the 300 GeV and the width below 1 GeV (as for tan (3 ~ 10) 
was assumed. For the Z/^* , the cms energy of the produced tt pair was 
taken in the range 300 it 10 GeV. A simple selection was applied. The 
minimal transverse momenta of the tt'.s were required to be above 15 GeV 
and the pseudorapidity \rj\ < 2.5 [9]. In Fig. 5 we can see the visible effect 
of spin-correlations. Similarily as in case of Fig. 3, the distribution for the 
Z/7* has more profounded shoulders than for the Higgs boson, due to spin 
correlations. 
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5. Summary 

We have discussed the algotithm for interfacing the r lepton decay pack- 
age TAUOLA with "any" production generator to include effects due to spin 
in the elementary Zj^* t^t^ process [4]. The invariant mass distribu- 
tions presented here, sensitive to the r+r^ spin correlations, can possibly 
be used for the MSSM Higgs boson searches at LHC to enhance sensitivity 
of the signal or to verify the hypothesis of the spin zero nature of the Higgs 
boson. This code is publicly availale from the address [12]. 
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